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Abstract 

In this report, generalized wave breaking equations are developed using 
three dimensional fully nonlinear extended Boussinesq equations to encompass 
rotational dynamics in wave breaking zone. The derivation for vorticity distri- 
butions are developed from Reynold based stress equations. 
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1 Introduction 

Wave breaking is one of the most complex phenomena that occurs in the nearshore 
region. During propagation of wave from deep to shallow water, the wave field is 
transformed due to shoaling and close to the shoreline, they become unstable and 
break. In the process of breaking, energy is redistributed from fairly organized wave 
motion to small scale turbulence, large scale currents and waves. 

It has been shown by numerous researchers that Boussinesq-type equations for 
varying water depth can describe nonlinear transformation in the shoaling region 
quite well. In the last couple of decades, a lot of research effort has gone into im- 
proving the predictive capability of these equations in the intermediate water-depth 
and close to the surf zone (see e.g. Nwogu [1983], MadsenjH] , Wei [1995]). It was 
established that to extend the validity of these equations to the deep water, higher 
order dispersive terms will have to be retained, and to improve the predictive capa- 
bility close to breaking, the nonlinear terms will all have to be retained. However, 
to model wave breaking, these models use additional terms that artificially added to 
the momentum equation, which would then reproduce the main characteristic of a 
breaking wave, i.e. the reduction in wave height. For example, wave breaking in FUN- 
WAVE (FUNWAVE is based on the model described by Nwogu [1993]) is modeled 
by introducing momentum mixing term developed by Kennedy et al [1999]. 

Most progress have been done for potential flow, starting with the work of Nwogu 
[1993] and Madsen [1983]. Some work have also been done to address partially ro- 
tational flows by Shen [2000] . In the breaking region and in the surf zone, the wave 

1 Center of Higher Learning, University of Southern Mississippi, Bldg 1103, Stennis Space Center, 
MS 39529. Email : rdutt@orca.st.usm.edu(corresponding author) 

2 GeoResources Institute, Mississippi State University, Bldg 1103, Stennis Space Center, MS 
39529. Email : veeramon@gri.msstate.edu 



1 



breaking introduces vorticity into the fluid. To address this problem, Veeramony & 
Svendsen [2000] derived breaking terms in Boussinesq equation assuming flow as two- 
dimensional rotational flow. Here, the breaking process is modeled by assuming that 
vorticity is generated in the roller region of the breaking wave and solving vorticity 
transport equation to determine the distribution of the vorticity. This naturally intro- 
duces additional terms in the momentum equation which causes wave height reduction 
as well as changes in the velocity field. However, since this model is based on stream 
function formulation, it cannot be trivially extended to three-dimensional flow. The 
phenomena of wave breaking in Boussinesq equations are being modeled using quite 
few techniques which can preserve the wave shape as well as include energy dissipa- 
tion mechanism. Shen [2000] developed a generalized form of Bousinesq equation in 
3D vortical flow field with arbitrary vorticity distribution up to 0(/z 2 ). But he did 
not describe momentum transport equation with full description of rotational flow. 
Recently, Zou et al [2004] addressed the problem by including the higher order terms 
in Boussinesq equation in 2D flow. This model solves to vorticity distribution based 
on the parametric form taken form surface roller data. In this paper, we try to de- 
velop a general form for breaking term for fully nonlinear set of Boussinesq equations 
for three dimensional vortical flow field near surf zone region. Derivation of breaking 
term from Reynold stress based vorticity transport equation was also developed to 
describe rotational field as a complete model of Veeramony [2000] . 

The paper is organized as follows: Section 2 discusses the basic governing equa- 
tions for continuity and momentum with boundary conditions. Section 3 describes 
the equation for horizontal and vertical velocity distribution for potential and rota- 
tional components. In section 4, the breaking term is derived for velocity transport 
equation for fully nonlinear case and solved vorticity transport equation analytically 
from fourier series expansion. In last section, results were discussed with conclusion. 

2 Basic Equations 

We consider a three-dimensional wave field with free surface r](x, y, t) propagating over 
a variable water depth h(x,y). As we are primarily concerned with wave breaking, 
we only consider here wave propagation in shallow water. Wave in this region can be 
characterized by two non-dimensional parameters 5 = a/h and /i — h/l where a is the 
characteristic wave amplitude and I the characteristic wave length. The parameter \x 
is a measure of frequency dispersion and 8 that of the nonlinearity of the wave. In this 
study, since we are only considering shallow water waves, we only have to consider 
weakly dispersive waves (upto 0(/z 2 )) but have to retain all nonlinear terms. 
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In this paper, the variables are non-dimensionalized using following scaling: 



x = x/l, y = y/l, z = z/h, t = t\/gh/l, (la) 
it = (^\fgh s J u, v — (^\fgh s J v, w — (5u\fgh^ w (lb) 

where the " represents the dimensional variables, g is the acceleration due to gravity, 
u and v are the horizontal components of the velocity in the x and y directions respec- 
tively, w is the vertical velocity. We start with the Eulerian equations of continuity 
and momentum in nondimensionalized form for velocity field u = (u,v,w) as: 

du . du _ du . du dp n 

— + 5u— + 5v— + 5w— + tt = 2 

at ox ay oz ox 

dv . dv . dv r dv dp „ 

— + Sv— + Sv— + 5w— + -f = 3 

at ox oy oz Oy 

5 ^ + 5 » U d^ + 5 » V dy- + 5 » W dz- + 5 d-z + 1 = ° (4) 

Since the fluid flow is rotational, we also have three dimensional vorticity field 
s = (s x ,s y , s z ) in the fluid defined as 

V x u = s (5) 
where y = (d/dx,d/dy,d/dz). The continuity equation then becomes, 

V» + ^ = (6) 

Here y • u — (du/dx,dv/dy). The above equations satisfy two boundary conditions 
for velocity at bottom and at free surface. At the free surface z = r)(x,y,t), since 
particles are free to move with fluid velocity, the kinematic boundary condition is 

w v = u v ■ yr] + — (7) 

and at bottom z = —h(x, y) 

w b = -u b ■ y/i (8) 

where = (u^jV^) is two component horizontal surface velocity. = (i] x ,i]y), 

Xjh = (h x ,h y ) refer to horizontal derivative with respect to x and y in all subse- 
quent calculations. The horizontal component for vorticity field s = (s y , — s x ) can be 
described as, 

du o , . 

— - a xjw = s (9) 
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with u = (u, v) as two component horizontal field whereas vertical component of 
vorticity expressed as 

_ s = du_dv_ , 1Q . 
2 dy dx 

This is straightforward calculation from equation (6) and (8) which is the begin- 
ning equation in three dimensional vorticity field formulation. 

d 2 iv 

pt\ 2 w + = ~ V • s = S w (11) 

w represents the vertical velocity of the flow. 

In the above equation, once w solved, horizontal component of velocity u, v can 
be solved from vorticity relation. In weakly hydrostatic case ( < /i 2 <C 1 ), solution 
is typically obtained from iterative perturbation procedure with successive correction 
term up to fi 2 . 

In case of breaking waves where vorticity is very strong, so (du/dz ~ 0(1)). We 
assume solution as, u = u Q + fi 2 ui + 0(/i 4 ) and w = w + fi 2 Wi + 0(/x 4 ) for horizontal 
and vertical velocity component. 

Under this assumption, Poisson equation becomes 

S2W ° (12) 



dz 2 



d 2 wi 
~dz^ 



d 2 w d 2 w 



dx 2 ' dy 2 J ^ 13 ^ 

wo, w\ can be calculated from bottom boundary conditions using equationn (7) sep- 
arately where the boundary conditions are, 



and 



w b0 = ~u b0 ■ \/h (14) 
Wbi + u bl ■ \/h = (15) 



at bottom boundary z = -h 

Since at any other depth z = z r , w is constrained by continuity equation only, so 
the equation follows 

-r^\z r = -X] -u m \ Zr + — -X]z r \ Zr (16) 

where u m is velocity at any arbitrary depth z r . . In Boussinesq type equation, one 
may take depth average or any intermediate velocity for horizontal velocity between 
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bottom and free surface as reference velocity. In the wave breaking zone where the 
vorticity is developed non uniformly, the equations become simpler with the choice of 
depth average velocity which includes contribution from surface vorticity gradient. We 
assume solution for velocity comes also from rotational contribution due to vorticity at 
the wave surface. So the velocity has both potential as well as rotational component, 
u = u p + u r , w = Wp + w r We solve w , W\ and u , u\ at any depth z r 



dw 
dz 1 



= -V • («™ - z r s\ Zr ) + z r y ■ s (17) 

r<9u>l. 



dz 



= [V^oW • X/z. 
du 



and 



dz 
dui 



(19) 

V^o (20) 



dz 

with boundary condition [uo] Zr = u r and [ui] Zr = Equations (4) - (16) form basic 
shallow water Boussinesq equations. 



3 Equation for horizontal velocity 

In the surf zone, vorticity grows very strongly as a non uniform function over depth. 
Following Shen [2000], we define reference velocity as u = u + Au — r/s^ in terms 
of depth average velocity u and magnitude of vorticity at free surface s v with the 
assumption of y • s ^ 0. we set here z r = i] as linear calibration for z r — r(rj + h) — h 
does not hold here in presence of nonuniform velocity as wave dispersion properties 
change both spatially and temporally with vorticity. And boundary condition can be 
set as 

dw . , . . . 

-fa\v = V ■ U + V(V ■ Sri) ( 21 ) 

Integrating equation (9) from bottom to surface and applying boundary condition 
to (16) we get wo as , 

w = w b0 - (-y -u + i]\/ ■ s v )H z - S w0 (22) 

where 

S w o= / (- V -s)dzdz 
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is the vertical velocity distribution generated by horizontal divergence of vorticity 
added to surface velocity. 

Now, once wq is calculated, u\ can be calculated from eqn (15) with surface 
boundary condition [uo] v = u rn and [ui] v = 

Finally, we calculate horizontal velocity as 



u 



sdz + fi 2 (S w i - S w i) 
+ ^(H 2 v -H z 2 )s7(vu-vV-s v ) 

+ /j, 2 (H v - H z ) [v((u + vs v ) -Vh) + (V -u-vV s v ) V h] 

+ 0(/i 4 ) (23) 

which on averaging over depth yields, 

u 2 

u = u v - Au + —H 2 v (V ' u ~ V V - s v) 

- Y H rilv(u + r]s v ) ■ sjh - (y -u - T) v Sr,) V h] +0(/i 4 ) (24) 

Au = jj- j_ h Au(z)dz is the average surface velocity contribution due to vorticity 
and it is significant for suspended sediment particles in the flow. The term Au(z) = 

sdz is the change due to depth variation of vorticity S. The total water depth H z 
and surface elevation H v are taken as H z = z + h and H v = r\ + h 

The contribution for velocity has and rotational component apart from potential 
due to vorticity generation. After we redefine H v = d and z = H z /H v , we express 
potential and rotational component up to order 0(/i 2 ) as 

fi 2 1 

u p (z) = u p + —(-- z 2 )d 2 v(V-«p) 

+ ^ 2 {\- z)d[\j{u p - \jh) + {xj-u p ) xjh] (25) 



u r (z) = u r - Au(z) + i]s v + fi 2 {S w i - S w i) 
~ y(^VlV(V-s,) 

- ^{\~ z)d7][(\j ■ s v ) xjh- v(s„- \jh)\ (26) 
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Similar expressions for vertical velocity are 



w p (z) = -(h + z) v -Up(z) 

u 2 1 

= - V {(h + z)u p ] - ^-(- - z 2 ) V \d 2 {h + z)(v(V ' «p)] 

- ^-2)[V'(Hz)(v(« P -V'i) + (V-« P )v'i)] (27) 

•uv(z) = - v '(h + z)u r - v ■ \{h + z)rjs v ] 

- ^v(H.)(^y V (ivg] 

+ V 2 V -[(\- z)d[v(vsr,- Vh)-(wsr,) Vh] (28) 

4 Breaking Model [fully nonlinear case] 

Conventional time dependent Boussinesq equations for surface wave height and con- 
sequent breaking term calculation are very straight forward and published previously 
in case of irrotational waves. Here we take up fully nonlinear calculation as vortic- 
ity becomes a large fraction of water depth in the surf zone or shoaling waves. So, 
while developing Boussinesq equations for horizontal momentum, we retain up to 
order 0(<5 2 ) and 0(<5/i 2 ) in our fully nonlinear calculation. Fully nonlinear Boussi- 
nesq equations for long wave have been derived by Mei [1983] for flat bottom and 
by Wei et al [1995] for variable bottom surface in case of irrotational wave. Shen 
[2000] addressed problems in developing generalized three dimensional irrotational 
propagating wave field to include rotational motion in general did not describe the 
vorticity breaking terms. For horizontal propagation of waves, the three dimensional 
problem can be reduced in terms of two horizontal velocity by integrating over depth 
and retaining up to order 0(5 2 ) and 0(5fi 2 ) As horizontal velocity is governed by 
momentum equation at the surface 77 by, 

Du. ,Dw, 

^l, = (^l, + l)Vr, (29) 

In the surf zone region of sloping beach, waves break due to high vorticity and the 
breaking of wave later being converted to turbulence. So horizontal variation of water 
depth h(x, y) must be considered in this case. We express surface propagation equa- 
tion in terms of average velocity description and total time derivative of horizontal 
momentum can be written as, 

Du du 
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where surface velocity is given by, 



u v = u + riS v - —d 2 v ( V • u ~ V V -s„) 



+ ~ Am I-^ + r 7 s r?) • X/h) - (V ' u ~ V V -Stj) V h] (31) 



i 2 
2 

We consider \jH v = \/rj + \/h for wavy bottom 



-Dm . (9m <9s„ _ _ u 2 , 9 r , (9m <9s„ . 

+ u-\7(\7-u-ri\7-s v ) + —d[\7(— + ri-^)-\7h 
du ds 

+ M- V ] + 0(u 4 ) (32) 

This long wave momentum equation upon simplification over flat bottom case 
can be compared to the one derived by Shen [2000] The vertical velocity can be 
obtained similarly, 

Dw . dw . dw . , noN 

-jxU = -^\v + u vVWr, + w—\ v (33) 

So, we can write the horizontal momentum equation as, 

du _ u 2 ,, r , . , du, 

— + u-Xju + Xjr 1 = —d 2 {xj{xj-u)- XJU + XJ{XJ-—) 

du du 
+ (« • V(V ' «)} - " 2 4V • - d 2 V (V • -Qj;) + (V • u v ) 2 

- m-v(V'«)}V^ (34) 

u is defined in previous section. In contrast to the result by Shen [2000], additional 
contribution factor here arises from vorticity variation which is significant for surf 
zone wave. Wei et al [1995] also breaking term for irrotational long wave momentum 
equation over a variable bottom wave. The intermediate depth velocity z a is being 
used there proportional to h instead of depth average velocity used here which may 
not be valid inside the fluid. The use of z r . in our approach avoids this difficulty 
Finally we try to generalize equation by solving vorticity from vorticity transport 
equation in next section. 
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5 Vorticity transport equation in breaking zone 



Madsen and Svendsen [1983] used a cubic vertical distribution of rotational velocity 
based on roller jump data which can not considered in three dimension case as it 
is not guaranteed to bring accuracy in the simulation. So we try to solve vorticity 
function from Reynold stress based equation. 

9u , . 1 , or . 

— + (u ■ \7)u = --VP (35) 

Taking the curl on both sides and use vorticity function s = y x u we get, 

ds 

— -(s- v)s + (u- y)s = z/y 2 s (36) 

(s ■ y)s is "vorticity stretching" factor due to change in gradient in vorticity. This 
term leads to change of rotation of material particles present in the flow to the beach. 
Contribution of this term can not be incorporated from two dimension roller jump 
data. 

We generalize the equation in three dimension as 

ds ds ds ds du dv dw 

»T + du TT + dv 7T + dw 7T ~ ds 7> ds 7> ^ s ^~ 

dt dx dy dz dx dy dz 

2 d 2 s 2 d 2 s d 2 s 
= U[fi d? +fI dy- 2+ ^ ] (37) 

After changing the variable from (x,y,z,t) to wave following coordinates (x, y, a, t), we 
write the vorticity equation as 

ds 5a drj ds 5 dw ds 

dt - iV ' S) ~ iV ' U) ~W^) [S ^~ W da- ] 



V W V 2 s + 77 — x r ^] + 0(/2 2 ) + 0(h x ) + 0(h y ) (38) 



(h + 5r]) yv "da (h + 5r]) yv "da 

1 d 2 s 

(h + Sri) da 2 

The boundary conditions in new coordinate system are, 

s(a = l,t) = s(x,y,t);s(a = 0,t) = 0; s(a, t = 0) = (39a) 

After we redefine s = Q + auj s , which transforms the equation to which is easier 
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to solve: 



dQ du s a di] dr]dVt 

~dt + °~dt~ JhT~6rf) ~dt ~ ~di~da 

. . „ x 5a . .dQ 5a 2 . . 

+ *" (v ' n) -(ftT^ (v ''' ) *-(ftT*rt (v ''' ) 

5u dfl 5w 
+ (h + 5ri)da~ + (h + 5r ] ) UJs 
5fl dw 5au s dw 



(h + 8rj) da (h + 5r]) da 

. . . . 5 2 afl du 5 2 afl . „,du 

8 2 au s , .du 

+ < 40) 

with new boundary, 

fl(a = l,t) = Ofl(a = 0, t) = (41) 

with initial condition fl(a, t = 0) = 0. This additional equation can be solved numer- 
ically as done by Briganti et al [2004] for the two dimensional case or an analytical 
solution can be formulated as shown by Veeramony & Svendsen [2000]. The analyt- 
ical solution can be calculated by assuming fl = + 5uj^ which gives first and 
second solution as 0(1) Problem 

duV> du s v d 2 uj^ , , 

where the solution is 

assuming —a^f- = ^2 n=1 °°F^ sinn7r<7 And to solve lo^\ assume u/ 1 ) = T, n Gn^ sinnna 
which gives zeroth order solution as 

G W = (-I)-— f ^±e n2 ^-t) dT (44) 
nrr J dr 

To consider 0(5) Problem 

djv v Q^m 



da h da 2 



= F (2) (45) 
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where 



adnduo^ a 2 dri a, .du;^ a 2 . . , m , . , 

F<> = ki^-i;i--^-^-^-n)+«(v^) (46) 



+ T- 



h da 

To solve above equation, assume = Y$ sinnna where solution becomes 



= 2 f F^e n2 ^< T - l) dT (47) 
Jo 



with 



F {2) = 2 [ F^sinnnada (48) 

The solution for vorticity s becomes, 

s = au) s + SiGn*- 1 -* sinmxcr + HiG n ^ sinnna (49) 

To solve breaking term, we need value of u s for boundary and eddy viscosity value 
as input data. 

6 Conclusion 

Finally we conclude here by developing a most generalized form of fully nonlin- 
ear Boussinesq equations for wave propagation in surf zone region with variable 
bathymetry with vorticity distribution from Vorticity Transport Equation(VTE). In 
this wave breaking zone, vorticity generated by the shear stress of current is very 
strong, so contribution to the surface velocity due to vorticity variation has signifi- 
cant contribution in fluid flow. These extra terms in generalized equation complicate 
the numerical technique as these terms are present in the equation in multiple form 
of equations for vorticity components which has to be solved in coupled solution tech- 
nique. Veeramony [2000] used simplied the formulation by taking constant eddy vis- 
cosity value but this oversimplified case may bring inaccuracy in calculation. Briganti 
et.al [2004] formulated a numerical technique scheme to solve VTE using generalized 
depth variable eddy viscosity v = u(x, y) in two dimension case. In three dimensional 
formulation, the nonlinear terms in the vorticity transport equation(VTE) will com- 
plicate the calculation and so proper numerical technique have to be developed. This 
work is under way. 
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